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ABSTRACT 

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton 
selected targets covering the 2 deg^ COSMOS field, with a yield of 485 high-confidence redshifts. 
The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS in- 
strument on the Magellan (Baade) telescope. We also include in the sample previously published 
Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We de- 
tail the observations and classification analyses. The survey is 90% complete to flux limits of 
/o.5-iokcV > 8 X 10"^^ erg cm~^ s~^ and i^g < 22, where over 90% of targets have high-confidence 
redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a 
description of the complete population to < 23. The corrected sample includes 57% broad emis- 
sion line (Type 1, unobscured) AGN at 0.13 < z < 4.26, 25% narrow emission line (Type 2, obscured) 
AGN at 0.07 < z < 1.29, and 18% absorption fine (host-dominated, obscured) AGN at < 2: < 1.22 
(excluding the stars that made up 4% of the X-ray targets) . We show that the survey's limits in X-ray 
and optical fiux include nearly all X-ray AGN (defined by Lo.s-iokcV > 3 X lO'^^ gj.g g-i) to z < 1, 
of both optically obscured and unobscured types. We find statistically significant evidence that the 
obscured to unobscured AGN ratio at z < 1 increases with redshift and decreases with luminosity. 
Subject headings: galaxies: active — galaxies: Seyfert — quasars — surveys — X-rays: galaxies 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) are the brightest 
persistent extragalactic sources in the sky across 
nearly all of the electromagnetic spectrum. Only in 
the relatively narrow range of infrared (IR) through 
ultraviolet (UV) wavelengths are AGN often outshone 
by stellar emission. Here the central engines can be 
dimmed by obscuring dust and gas while starlight, 
either direct or absorbed and re-emitted by dust, peaks. 
Historically, the largest AGN surveys have been based 
on optica l selection (e.g. [B OS, Schmidt & Grcenj 

1995; 
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ISDSS. Schneider et al.ll2007t ). Yet in both the local and 
distant universe, obscured AGN are generally thought 
to outn umber their u nobscured coun t erparts (e.g. 
Ivlaiolino I FRieki [1 995': Gi Ui et al.l [200l iSteffen et all 
|2004: Barger et al.l [2005; .Martinez-Sansigre et all 120051: 
lOaddi et all 120071: iTreister. Krolik fc Dullemo nd 2008), 
indicating that optical surveys probably miss the 
majority of AGN. A more complete census of AGN 
must use their X-ray, mid-infrared, and radio emission, 
where obscuration and host contamination are mini- 
mized. X-ray and mid-IR selected surveys do in fact 
reveal a far greater space density of AGN than optical 
selection: for example, the Chandra deep fields reveal 
AGN sky densities 10-20 times higher than those of 
optically sele cted surveys to the same limiting optical 
magnitudes (iBauer et all 12004 iRisaliti fc Ehdi (20041 : 
iBrandt fc Hasingeil I2005D . However, most X-ray and 
mid-IR surveys either have significantly smaller areas 
and numbers of AGN or are wide-area b ut sub stantially 
shallower than optical s urveys (e.g. ISchwope et all 
120001 : iLonsdale et al.l I2003D . Here we present a deep 
spectroscopic survey of AGN both without the biases of 
optical selection and over a relatively large field. 

The Cosmic Evolution Survey (COSMOS, 
iScoville et all [200l ^° is built upon an HST Trea- 
sury project to fully image a 2 deg^ equatorial field. 
The 590 orbits of HST ACS i-band observations have 
been supplemented by observations at wavelengths from 
radio to X-ray, including deep VLA, Spitzer, CFHT, 
Subaru (6 broad bands and 14 narrow bands), GALEX, 
XMM-Newton, and Chandra data. Here we present a 
complete spectroscopic survey of XMM-selected AGN 

1" The COSMOS website is |http: / /cosmos. astro.caltech.edu/ 1 
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in the COSMOS field. Most (601) target s have spectra 
taken with the IMACS spectrograph (jBigelow et al.l 
|1998[) on the MageUan te lescope, includ ing 282 spectra 
previously pubhshed by (|Trunip et al.l [2007). An ad- 
ditional 76 X-ray targets were excluded from IMACS 
observations because they already had SDSS spec- 
tra. For 134 of the targets with IMACS coverage, 
we additionall y acquired spec t ra wi th the Hectospec 
spectrograph (|Fabricant et all 120051 ) on the MMT 
telescope as ancillary data with extended blue coverage. 
In total, we were able to target 52% (677/1310) of 
the available z^g < 23.5 X-ray targets, resulting in 
485 high-confidence redshifts. The relevant observing 
strategics and configurations are described in detail in 
§2. We were 90% in assigning high-confidence redshifts 
to all spectral types at i^g < 22, with decreasing confi- 
dence, dependent on both redshift and spectral type, at 
fainter magnitudes. The IMACS spectroscopy campaign 
additionally targeted AGN candidates selected by their 
radio (VLA, 605 targets) and IR (Spitzer/IRAC, 236 
targets) emission, but these objects are not included in 
this study and will be presented in future work. 

We place this work in the context of other large X-ray 
AGN surveys in Figure [1] where the left panel compares 
the X-ray depth, areal coverage, and number of sources 
for various X-ray AGN surveys. The right panel of Fig- 
ure [T] s hows our flux limits wi th the customary "AGN 
locus" ()Maccacaro et al.l 119881 ). The depth of XMM- 
Newton in COSM OS most closely resembles the AEGIS 
([Davis et al.ll2007[) survey, with roughly the same number 
of X-ray targets in both despite their slight differences in 
area and X-ray depth. There exists no purely optical sur- 
vey to the depth of our spectroscopy (i^g < 23.5) with 
this number of spectroscopic redshifts. The AGN spec- 
troscopic campaign presented here is significantly deeper 
than large optical surveys like the 2dF Quasar Redshift 
Survey i2dF, Groom ct al. 2001 ) and th e Sloan Digital 
Sky Survey (,SDSS. Schneider et alll2007t l. In particular, 
we present targets ~60 times fainter than the main SDSS 
spectroscopy {g < 19.1), and ~20 times fainter than the 
deepest SDSS spectroscopy {g < 20.2) for quasars, and 
our spectroscopy reaches a (arbitrary) quasar/Seyfert 
boundary of = -23 at z ^3. Surveys like the VI MOS 
Very Deep Survey (VVDS, iGavignaud et al.ll2006[ ) may 
reach similarly faint magnitudes (z < 24 in VVDS) but 
have far fewer AGN (130 in VVDS). We additionally note 
that the Magellan AGN sample will eventually be aug- 
mented by ~300 X-ray AGN from the faint zCOSMOS 
surve y of galaxy redshifts with VLT/ VIMOS ()Lillv et all 
[200l . 

We discuss the analysis of the spectra in §3, includ- 
ing the methods for classifying the AGN and determin- 
ing redshifts. In §4 we characterize the completeness of 
the survey and discuss the populations of different AGN 
types. We use the sample to understand the X-ray AGN 
population in §5, and we discuss future projects using 
this dataset in §6. We adopt a cosmol ogy consistent 
with WMAP results (|Spergel et al.l[200l oi h = 0.70, 
Qm = 0.3, r^A = 0.7. 

Throughout the paper we use "unobscured" to describe 
Type 1 AGN with broad emission lines and "obscured" 
to describe X-ray AGN where the host galaxy light dom- 
inates the optical continuum. Thus we use "obscured 



AGN" to describe both spectroscopically-defined Type 
2 AGN (with narrow emission lines, classified as "nl" 
or "nla" in the catalog) and XBONGs (X-ray bright, 
optically norm al galaxies, classified as "a" in the cat- 
alog, see also [ConmstaTeTaD 120021 : iRigbv et al.l 120061 : 
ICivano et al.l 12007) . It is important to note that our 
designation as "obscured" does not necessarily describe 
the physical reason for the faint optical nuclear emis- 
sion: the AGN might simply be under- luminous in the 
optical instead of being hidden by obscuring material. 
Indeed, ma ny Type 2 AGN app ear to be unobscured in 
the X-rays (|Szokolv et al.ll2004l ). while broad absorption 
line (BAL) Type 1 AGN are typically X-ray obs cured 
(|Brandt. Laor. fc Will£|[2000t iGallagher et al.li2006l) . We 
also note that even our "obscured" AGN types have mod- 
erate X-ray luminosity and we are not sensitive to heav- 
ily X-ray obscured (e.g., Compton-thick, Nh ^ 1 x lO^'' 
cm~2) AGN which are too faint for our XMM-Newton 
observations. 

2. OBSERVATIONS 
2.1. XMM 

The COSMOS field has been observed with XMM- 
Newton for a total of 1.55 Ms at th e homogeneous 
vignetting-corrected de pth of ~ 50 ks ([Hasinger et al.l 
[2007: CappcUuti ct al. 2007 ll2008D . The final catalog in- 
cludes 1887 point-like sources detected in at least one 
of the soft (0.5-2 keV), hard (2-10 keV) or ultra-hard 
(5-10 keV) bands down to limiting fluxes of 5 x 10^^®, 
3.3 X 10~^^, and 5 x 10~^^ erg cm~^ s~^, respectively 

see ICappelluti eran[2007L [20081 . for more details). The 
detection threshold corresponds to a probability < 4.5 x 
10"^ that a source is instead a background fluctuation. 
The XMM fluxes have been computed converting the 
count-rate into flux assuming a spectral index F = 2.0 
and Galactic column density Nh = 2.5 x 10^° cm^ for 
0.5-2 keV and F — 1.7 and Galactic colu mn density 
Nh = 2.5 X 10^0 cm^ for 2-10 keV. Following iBrusa et all 
([2008f ). we exclude 24 sources which are a blend of two 
Chandra sources and 26 faint XMM sources coincident 
with diffuse emission (iFino guenov et al. 2008). We im- 
pose a brighter flux limit than the full catalog because the 
XMM-Newton observations were not complete until the 
3rd season (2007) of spectroscopic observing. Figure [2] 
shows the X-ray sensitivity for each of the three seasons 
of IMACS, revealing that the first two seasons (2005- 
2006) suffer from slightly shallower X-ray catalogs. The 
sample we use is limited to flux limits of the 50% XMM 
coverage area, which has only 186 few sources than from 
the limits of the entire XMM coverage. The sample in- 
cludes 1651 X-ray sources detected at fluxes larger than 
1 X IQ-^^ cgs, 6x 10-15 cgs, 1 X lO-i'^ cgs, in the 0.5-2 keV, 
2-10 keV or 5 -10 ke V bands, respectively, as presented by 
iBrusa et al.l (I2008D. 

iBrusa et alllM ) associated the X-ray point sources 
with optical counterparts using the likelihood ratio 
technique to match to the optical, near-infrared (K- 
band) and mid-in frared (IRAC) photometric catalogs 
([Capak et al.l[2007h . The images for the XMM-COSMOS 
subsample additionally covered by Chandra observa- 
tions were matched to the Chandra/ACIS images by 
visual insp ection (Elvis et al] 1200 8: Pucc etti et al.ll2008l : 
ICivano et al . 2008). We use the COSMOS Chandra ob- 
servations for reliability checks only, since it covers only 
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the central 0.8 deg^ and is still undergoing basic analyses. 

Of the 1651 sources in the XMM-COSMOS catalog 
described above, 1465 sources have an unique/secure op- 
tical counterpart from the multiwavelength analysis with 
a probability of misidentification of < 1%. For an addi- 
tional 175 sources, there is a second optical source with 
a comparable probability to be the correct counterpart. 
Because the alternate counterpart shows comparable op- 
tical t o IR properties (and comparable photometric red- 
shifts, ISalvato et all l2008f ) to the primary counterpart, 
the primary counterpart can be considered statistically 
representative of the true counterpart for these 175 X- 
ray sources, and we include the primary counterparts in 
the target sample. Eleven sources (outside the Chandra 
area) remain unidentified because they had no optical or 
infrared counterparts (i.e., their optical/infrared counter- 
parts were fainter than our photometry) . We designated 
the 1310 optical counterparts with ijg < 23.5 (from the 
CFHT) as the X-ray selected targets for the spectroscopic 
survey. 

2.2. Magellan/IMA CS 

The bulk of the spectroscopic data comes from obser- 
vations with the Inamori Ma gellan Areal C amera and 
Spectrograph (IMACS, Bigelo w et all 119981 ) on the 6.5 
m Magellan/Baade telescope. The IMACS field of view 
is 22'30" X 21'10" (with only 10% vignetting at the ex- 
treme chip edge), requiring 16 tiled pointings to fully ob- 
serve the entire 2 deg^ COSMOS field as shown in Figure 
[3l We observed these 16 pointings over the course of 26 
nights (18 clear) through three years, as detailed in Table 
[TJ The total exposure time for each pointing is 4-6 hours 
(shown in Table [T] and Figure [3]). Henceforth we refer 
to each pointing by its number in Table [T] and Figure [3l 
We were able to simultaneously observe 200-400 spectra 
per mask: generally ~40 of these were the X-ray targets 
described here (shown in the last column of Table [1]), and 
the additional slits were ancillary targets to be described 
in future work. We were generally able to target ~50% 
of the available i^g < 23.5 X-ray targets in each tiled 
IMACS field, or 601/1310 X-ray targets over the 2 deg^. 

All IMACS spectra were obtained over the wavelength 
range of 5600-9200 A, with the Moon below the hori- 
zon and a mean airmass of 1.3. We used the 200 1/mm 
grism in the first year and a 150 1/mm grism designed 
and constructed for COSMOS in the second and third 
years. The lower-resolution 150 1/mm grism had a reso- 
lution element of 10 A. Since all observed broad line AGN 
had line widths > 1500 km and all observed narrow 
line AGN had line widths < 1000 km s~^, the resolution 
of the grism was sufficient to distinguish broad and nar- 
row fine AGN. The gain in S/N from 200 1/mm to 150 
1/mm was only marginal, but the 150 1/mm grism al- 
lowed for a maximum of 400 slits per mask, ~35% more 
than the maximum 300 slits per mask for the 200 1/mm 
grism. The slits were 11" x 1" (55 x 5 pixels), though 
only 5'.'4x 1" of the slit was cut, so that an extra adjacent 
5'.'6 was reserved as an "uncut region" to accommodate 
"nod-and-shuffle" observing (see below). We attempted 
to observe each mask for 5 or more hours, which achieves 
high completeness of AGN rcdshifts at z^g ~ 23, al- 
though as Figure [3] shows this was not always achieved. 
We estimate the impact of the nonuniform spectroscopic 



depth on the sample's completeness in §4.1. 

We observed using the "nod-and-shuffle" technique, 
which allowed for sky subtraction and fringe re- 
moval in the red up to an order of magnitude 
more precisely than conventional methods. The gen- 
eral principles of nod-and-sh uf fle a re described by 
iGlazebrook fc Bland-HawthornI ()200lD . and our ap- 
proach is detailed in Appendix 1 of [Abraham et al.l 
(2004). Briefly, we began observing with the target ob- 
jects offset from the vertical center of the cut region, 1/3 
of the way from the bottom to the top (that is, 1'.'8 from 
the bottom slit edge, and 3'.'6 from the top edge of the 
cut region and the cut/uncut boundary). After 60 sec- 
onds we closed the shutter, nodded the telescope by I'.'S 
(9 pixels) along the slit, and shuffled the charge to the 
reserved uncut region. The object was then observed for 
60 seconds in the new position, 2 /3 of the way from the 
bottom to the top of the cut region (3'.'6 from the bottom 
and 1'.'8 from the top). We then closed the shutter, nod- 
ded back to the original position, and shuffled the charge 
back onto the cut region on the mask. This cycle was 
repeated (typically 15-20 times) with the net result that 
the sky and object had been observed for equal amounts 
of time on identical pixels on the CCD. Nod-and-shuffle 
worked well while the seeing was ^ 1", which was true 
for all observations. 

To extract and sky-subtract individual 2D lin- 
ear IMACS spectra, we used the publicly avail- 
able Carnegie Observatories System for MultiOb- 
ject Spectroscopy (with coincidentally the acronym 
"COSMOS," written by A. Oemler, K. Clardy, D. 
Kelson, and G. Walth and publicly available at 
http://www.ociw.edu/Code/cosmos). We combined the 
two nod positions in the nod-and-shuffle data, then co- 
added the individual 2D exposures of each pointing while 
rejecting cosmic rays as 4.5cr outliers from the mean of 
the individual exposures. Wavelength calibration was 
performed using an He/Ne/Ar arc lamp exposure in each 
slit. The 2D spectra were extracted to ID flux-calibrated 
spectra using our own IDL softw are, adapted from the 
ispec2d package (Mou stakas fc K cnnicutt 2006). While 
flux calibration used only a single standard star at the 
center of the IMACS detector, we estimate by eye that 
vignetting has < 10% effect on the spectral shape or 
throughput across the field, in agreement with the pre- 
dictions of the IMACS manual. 

IMACS spectra can be contaminated or compromised 
in several ways, including 0th and 2nd order lines from 
other spectra, bad pixels and columns, chip gaps, poorly 
machined slits, and cosmic rays missed during co-adding. 
To eliminate these artifacts, we generated bad pixel 
masks for all ID spectra by visual inspection of the cal- 
ibrated ID and 2D data. The nod-and-shuffle 2D data 
were especially useful for artifact rejection: any feature 
appearing in only one of the two nod positions is clearly 
an artifact. Pixels designated as bad in the mask were 
ignored in all subsequent analyses. 

We show 10 examples of IMACS spectra in Fig- 
ures 13] and El These spectra are representative of 
the targets in the survey. Each of these spectra are 
smoothed by the 5-pixel resolution element. We dis- 
cuss each object below, with the spectral classifica- 
tion, confidences, and redshift algorithms detailed in 
§3. Briefly, Zconf = 3,4 refer to high confidence and 
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■Zconf = 1,2 are lower confidence guesses (but see also §3.1 
for the subtleties in confidence assignment). All spec- 
tra are publicly available on the COSMOS IRSA server 
(|http://irsa.ipac.caltech.edu/data/COSMOS/). 

1. COSMOS J095909. 53+021916. 5, i^^^ = 20.05, z = 
0.38, Zconf — 4: This is a low redshift Type 1 
Seyfert. The emission lines are bright and easily 
identified. 

2. COSMOS J095752. 17+015120.1, 4b = 21.00, z = 
4.17, Zconf = 4: This is a high redshift Type 1 
quasar. Lya is especially prominent along other 
broad emission features, and so this redshift is very 
reliable. 

3. COSMOS J095836.69+022049.0, ij^ = 23.04, z = 
1.19, Zconf = 4: We classify this target as a hybrid 
"bnl" object with both broad and narrow emission 
lines. The narrow [On] line is evident above the 
noise and strong broad Mgiiis also present. 

4. COSMOS J095756.77+024840.9, = 19.60, z = 
1.61, Zconf = 3: In this spectrum, a broad emission 
line is cut off by a detector chip gap. Identifying 
the broad feature as Mgll and the minor narrow 
emission line at ~6375A as [Niv] yields a good 
redshift, but we assign only Zconf = 3 because of 
the uncertainty from the chip gap position. 

5. COSMOS J100113. 83+014000.9, i^^ = 20.49, z = 
1.56, Zconf = 2: The blue end of this spectrum lies 
on a chip gap, and much of the red end is corrupted 
by second order features from another bright spec- 
trum on the mask. Only one broad emission line is 
present, and so while the target is clearly a Type 1 
AGN, the line could be either CIII] or Mgll. The 
redshift solution is degenerate and we assign only 

•^conf — 2. 

6. COSMOS J095821. 38+013322.8, ijg = 19.16, z = 
0.44, Zconf — 4: This spectrum contains several 
bright emission lines, and is clearly identified as a 
"nl" class object. This object has Lq.s-io kcv < 
3 X 10''2 and -2 < log fx/ fo < -1, and it is prob- 
ably a starburst galaxy (see §3.2 for our distinction 
between AGN and starbursts). 

7. COSMOS J095855. 26+022713. 7, = 22.07, z = 
1.13, Zconf = 4: This narrow emission line spectrum 
is faint, but the [On] emission feature has a strong 
signal above the noisy continuum. We assign this 
spectrum Zconf = 4 because there is no other plau- 
sible redshift solution for a single bright narrow 
emission line. The 2D spectrum (not shown) also 
reveals the emission feature in both nodded posi- 
tions, confirming that it is not a noise spike. This 
object is a Type 2 AGN with both Lq 5-10 kcV > 
3 X 10^2 erg s'^ and -1 < log fx/ fo < 1- 

8. COSMOS J095806. 24+020113.8, i^^^ ^ 21.26, z = 
0.62, Zconf = 4: We identify this spectrum as a 
hybrid "nla" object, since it has both narrow emis- 
sion lines and the absorption lines of an early-type 
galaxy. H/? is present only in absorption, and while 



half of the H+K doublet is on a masked-out region, 
the other line is present. 

9. COSMOS J095906.97+021357.8, = 21.11, z = 
0.76, Zconf = 4: This spectrum exhibits only ab- 
sorption lines and is classified as an early-type 
galaxy. The continuum shape and H+K doublet 
make assigning redshifts to these targets straight- 
forward. This object meets both of the X-ray emis- 
sion criteria of §3.2 and is an optically obscured 
AGN. 

10. COSMOS J095743. 85+022239.1, = 22.20, z = 
1.02, Zconf = 1: This spectrum is quite noisy. The 
single narrow line may be [0 11] , but it is not strong 
enough above the noise to reliably classify. Because 
its entire identification may be a result of noise, we 
designate this target as Zconf = 1- 

2.3. MMT/Hectospec 

We also obtained ancillary spect roscopic data using the 
Hcctospec fiber- fed spectrograph (jFabricant et al.l 120051 ) 
on the 6.5 m MMT telescope. The field of view for Hec- 
tospec is a 1 deg diameter circle, and in March 2007 
the COSMOS field was observed with two pointings of 
3 hours each, as shown in Figure [3l These pointings 
contained a total of 134 targets to zjg < 23.5 in 2.5" 
fibers. We observed with the 270 1/mm grism over a 
wavelength coverage of 3800-9200A, resulting m a res- 
olution of 3A. Because Hectospec is fiber-fed and the 
MMT has a brighter sky and generally poorer seeing than 
Magellan, MMT/Hectospec cannot reach targets as faint 
as those reached by Magellan/IMACS. Therefore we use 
Hectospec observations only as ancillary data on targets 
which already have IMACS spectra. 

The MMT/Hectospec observations were designed pri- 
marily to double-check the redshifts derived from IMACS 
spectra by adding the bluer 3800-5600A wavelength 
band. Figure [6] shows the observed peak wavelength with 
redshift for the strong broad emission line in Type 1 
AGN. With IMACS, the hmited red wavelength range 
means that broad line AGN at 0.4 < z < 1.9 and 
2.3 < z < 2.9 will have only one observed broad line, 
as shaded in the figure. These potentially ambiguous 
redshifts can be resolved using the Hectospec spectra. 
Even for targets with non-ambiguous redshifts, the ex- 
tended wavelength coverage allows for consistency checks 
and additional line measurements. 

In Figure [7] we show two objects where a high- 
confidence redshift could be assigned only after Hec- 
tospec spectra were additionally taken. The first of 
these, 095801.45+014832.9, was assigned Zconf = 2 and 
an incorrect redshift of 1.3 before the Hectospec data al- 
lowed us to correctly resolve the degeneracy and assign 
Zconf = 4. The second object, 100149.00+024821.8, had 
been assigned the correct redshift from its IMACS spec- 
trum but only Zconf = 2, and the Hectospec data con- 
firmed the otherwise uncertain solution and allowed us 
to assign Zconf = 4. In general, the additional Hectospec 
spectra revealed that we were ~75% accurate in assign- 
ing redshifts to IMACS spectra with degenerate redshift 
solutions. (We were better than the 50% chance proba- 
bility because we were occasionally able to fit to minor 
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features, e.g. Fell/III complexes, weak narrow lines like 
[On] and [Niv], or general continuum shape.) 

We reduced the Hectospec data into ID Hnear spec- 
tra with sky subtraction, flux calibration, and cosmic 
ray rejection using the publicly available HSRED soft- 
ware (written by R. Cool). We also used HSRED to 
apply an artificial flux calibration to correct the spec- 
tral shape, then flux calibrated the spectra using a mean 
correction from objects with both Hectospec and IMACS 
data. From the resultant spectral shape we estimate that 
this technique has flux errors in the blue and red ends of 
the spectra as large as '^20%. Since the analyses are lim- 
ited to flnding redshifts and performing simple line width 
measurements, errors of this magnitude are acceptable. 

2.4. SDSS 

We include 76 XMM X-ray sources with spectra pre- 
viously tak en as part of th e Sloan Digital Sky Sur- 
vey (SDSS, lYorket al.l[20nl . With redshifts already 
known, these targets were excluded from the main 
IMACS survey. These objects were selected using 
the publicly available SDSS Catalog Archive Server 
(http://cas.sdss.org/astro/) , which uses th e SDSS Data 
Release 6 (Adelman-McCa rthv et al.ll2008l ). Their wave- 
length coverage is 3800-9200A and they have a resolution 
of 3A. AU the SDSS targets are uniformly bright, with 
*AB ~ 21, and so they would certainly have been success- 
fully observed with IMACS had their redshifts not been 
previously known. Including these SDSS targets does not 
introduce any new incompleteness or complication to the 
sample. 

3. SPECTRAL ANALYSIS 

Our program as described above is largely motivated 
as an AGN redshift survey. We especially seek Type 1, 
Type 2, and optically-obscured (host-dominated AGN), 
though we also find a small contaminant fraction of local 
stars and star-forming galaxies. We attempt to separate 
the population of obscured AGN from star-forming and 
quiescent galaxies using X-r ay and optical colo r diagnos- 
tics. A companion paper (!Trump et al.l I2008D presents 
basic line measurements and estimates of black hole mass 
for the Type 1 AGN. 

3.1. AGN Classification 

We used three compos ite spectra from the Sloan Dig- 
ital Sky Survey (SDSS, lYork et al.l |2000D as templates 
for classifying the objects and determining their red- 
shifts: a Type 1 (b road emission line) AGN composite 
from 2204 sources (Vanden Be rk et al.l[200ll ). a Type 2 

i narrow emission line) AGN composite from 291 sources 
Zakamska et al.ll2003)) . and a red g alaxv composite from 
965 sources ( Eisenstein et al.ll2001[ ). The three template 
spectra are shown in Figure El We found that the Type 
2 AGN composite gave accurate redshifts for both star- 
forming galaxies and AGN with narrow emission lines. 
The red galaxy template was likewise accurate for a va- 
riety of absorption line galaxies, ranging from old stel- 
lar systems with strong 4000A breaks to post-starburst 
galaxies. Objects showing a mixture of narrow emission 
lines and red galaxy continuum shape and absorption fea- 
tures were classified as hybrid objects. We did not use a 
particular template for local stars, but stars ranging in 



temperature from 0/B to M types were easily visually 
identified. 

To calculate redshifts we used a cross-correlation red- 
shift IDL algorithm in the publicly available idlspec2d 
package written by D. Schlegel^^. This algorithm used 
a visually-chosen template to find a best-fit redshift 
and its associated la error. As discussed in §2.2, all 
masked-out regions were ignored in the redshift deter- 
mination. Note that the redshift error returned is prob- 
ably underestimated for objects with lines shifted from 
the rest frame with respect to each other, as is often 
the case between high-ionization (e.g. Civ) and low- 
ioniza tion (e.g. Mgli) bro ad emission lines in Type 1 
AGN (jSulentic et al.l 120001 ). We manually assigned red- 
shift errors for 6% (41/677) of objects where the cross- 
correlation algorithm failed but we were able to visually 
assign a best-fit redshift. 

Each object was assigned a redshift confidence accord- 
ing to the ability of the redshifted template to fit the 
emission lines, absorption lines, and continuum of the ob- 
ject spectrum. If at least two emission or absorption lines 
were fit well, or if at least one line and the minor con- 
tinuum features were fit unambiguously, the redshift was 
considered at least 90% confident and assigned Zconf = 4 



(64% of objects). Objects of Zconf = 3 (8% of objects) 
have only one strong line feature with a continuum or 
second less certain feature that make their assigned red- 
shift likely but not as assured. We assign Zconf = 2 (8% 
of objects) when the spectrum exhibits only one broad or 
narrow feature and the calculated redshift is degenerate 
with another solution. Objects of Zconf — 1 (5% of ob- 
jects) are little more than guesses, where a sole feature 
is present but has little signal over the noise, such that 
even the spectral type classification is uncertain. (It is 
notable, however, that the nod and shuffle observations 
helped to resolve real features from noise, since real fea- 
tures must occupy both nodded positions on the CCD.) 
If the signal-to-noise of the object spectrum was too low 
for even a guess at the redshift or spectral type, it was 
assigned Zconf = (13% of objects). We additionally as- 
sign Zconf = —1 to 13 targets with "broken" slits, severely 
contaminated by second order lines or mask cutting er- 
rors. In total, we were unable to assign redshifts for 15% 
of targets. Duplicate observations with IMACS and Hec- 
tospec indicate that redshift confldences of 4, 3, 2, and 
1 correspond to correct redshift likelihoods of 97%, 90%, 
75%, and 33%, respectively. These duplicate observa- 
tions are mainly estimated for brighter targets, however, 
and so the true likelihoods may be slightly lower. In to- 
tal, we were able to classify 573 spectra with Zconf > 
and we designate the 485 spectra with Zconf = 3,4 as 
"high-confidence" objects. We discuss individual targets 
spanning the classification types and confidence levels in 
§2.2 and shown in Figures |4] and El 

All of the objects observed in the sample are presented 
in Table [2l The classifications are as follows: "bl" for 
broad emission line objects (Type 1 AGN), "bnl" for 
objects with both broad and narrow emission (possibly 
Type 1.5-1.9 AGN), "nl" for narrow emission line ob- 
jects (Type 2 AGN and star- forming galaxies), "a" for 
absorption line galaxies, "nla" for narrow emission and 
absorption line galaxy hybrids, and "star" for stars (of 

publicly available at jhttp:/ /spectro.princeton.edu/idlspec2d-install.html| 
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varied spectral type). We further classify narrow emis- 
sion and absorption line spectra as AGN or inactive in 
§3.2 below. In total, 50% (288/573) of the classified tar- 
gets were designated "bl" or "bnl," 30% (171/573) were 
"nl" or "nla," 16% (92/573) were "a," and the remaining 
4% (22/573) were stars. Objects with a question mark 
under "Type" in Table [5] have too low signal-to- noise to 
venture a classification, although many of these objects 
are unlikely to be Type 1 or 2 AGN for reasons we dis- 
cuss in §4. As mentioned above, objects with Zconf = 1 
may be incorrectly classified. 

We summarize the efficiencies, from X-ray sources to 
targeting to redshifts, in Table [3l 

3.2. AGN, Starhursts, and Quiescent Galaxies 

We use the following X-ray emission diagnostics to 
classify the 485 high-confidence extragalactic objects as 
AGN: 

io.5-10 kov > 3 X lO''^ erg s'^ (1) 

-1 < log/x//o < 1, where log fx /fo = log(/x)+«AB/2 

(2) 

Each of these criteria have been shown by several 
autho rs to reliably (albeit conservativel y) select AGN 
(e.g.. [Hornschemeier et aLir2 001: Alexa nder et al.l|2001l : 
iBauer et al.ll2004HBundv et al. 2007), although it is im- 
portant to note that bona fide AGN (e.g., LINERs 
and other low-luminosity AGN) can be much less X-ray 
bright than these criteria. Equation 1 is derived from 
the fact that purely star-forming galaxies in the local 
univer se do not exceed Lp 5_inkeV 9^ 3 x lO'*^ erg s~^ 
(e.g., iFabbianolfTosl iColbert et all 1200? ). The X-ray 
luminosities of the sources are shown in Figure [9] along 
with the X-ray flux limit. Nearly all of the Type 1 
AGN (marked as crosses) lie above the AGN luminos- 
ity threshold. Equation 2 is the tradit ional "AGN lo- 
cus" defined bv iMaccacaro et ahl (|1988f ). shown for the 
sample in Figure [TUl Objects marked with x's have 
-^0.5-iokcV > 3 X 10^^ erg s~^, revealing that the two 
methods heavily overlap, with 94% (405/432) of the ob- 
jects that satisfy one of the criteria additionally meeting 
both. Only 53 "nl" and "a" objects do not meet either of 
the X-ray criteria, leaving us with 432 X-ray AGN that 
meet either Eq. 1 or Eq. 2 and have high-confidence 
redshifts. 

Using either Equations 1 or 2 selects all of the spec- 
troscopically identified Type 1 AGN, but it still may 
exclude some obscured AGN. The source classifica- 
tion diagnostic diagrams, based on the optical emission 
line measurements Baldwin. Ph ilips, fc TerlevichI ()1981l . 
BPT), are usually quite effective in classifying narrow 
emission line spectra as star forming gal axies or Type 
2 AG N, with a sound theoretical ba sis (iKewlev et al. 
2001^ and use in m any surveys (e.g., 'Kauff mann et all 
2003; Trc monti et al. 2004). The BPT diagnostic uses 
ratios of nebular emission lines ([Oiii]A5007/H/3 and 
[Nii]A6583/Ha) to distinguish between thermal emis- 
sion from star formation and nonthermal AGN emission. 
However, there are two limitations to the BPT diagnostic 
that make it inapplicable to our sample. First, most of 
the object do not have the appropriate lines in their ob- 
served wavelength range: most of the "nl" objects are at 
higher redshift and we are limited by the spectral range of 



IMACS. In addition, accurately measuring the line ratios 
requires correcting for absorption in Ha and H/3 from old 
stellar populations. Because the spectra have low reso- 
lution and limited wavelength range, we are unable to 
accurately fit and account for stell ar absorption . 

The color-based diagnostic of iSmolcic et al.l (|2008af l 
can be used to further classify the narrow emission and 
absorption line spectra which do not satisfy Equations 
1 and 2 but are nonetheless AGN. This selection tech- 
nique is based on the a tight correlation in the local uni- 
verse between the emission line flux ratios utilized for 
the spectroscopic B PT selection and the galaxies' rest- 
frame optical colors (jSmolcic et al.l[20"08a[ ). The method 
h as been well cahbrated on the local SDSS/NVSS sample 
in ISmolcic et al .' (2008a) and sviccessfully applied to the 
radio VLA-COSMOS data (Smolcic ct al. 2008b). Fol- 
lowing [SmolliZI^lal] (|2008b[ ) the rest-frame color for the 
narrow line AGN was compute d by fitting each ga laxy's 
observed optical-to-NIR SED (jCapak et all l2007l ). de- 
redshifted using its spectroscopic redshift, with a library 
of 10 0,000 model s pectra (Bru zual fc Chario t 2003). 
5+5SKi8l cic et al.l ()2008aD show that the color diagnostic is a 
good statistical measure, but may not be accurate for in- 
dividual objects. So while it further indicates that 17/53 
objects are obscured AGN outside the X-ray criteria, we 
do not include these objects as AGN and note only that 
the sample of 432 high-confidence X-ray AGN as defined 
by Equations 1 and 2 probably misses at least ~17 ad- 
ditional objects. 

4. COMPLETENESS 

Our targeting was based solely on the available X- 
ray data and the optical flux constraint of i^g < 23.5. 
While we were only able to target 52% (677/1310) of 
the available X-ray sources, the spectra obtained were 
constrained only by slit placement and so represent a 
random subset of the total X-ray population. There- 
fore our completeness limits can be determined from the 
success rate for the spectroscopy, which is dependent on 
optical magnitude, object type, and redshift. We char- 
acterize and justify the flux limits in §4.1 and 4.2, as 
well as the more detailed redshift completeness in §4.3. 
Our goal is a purely X-ray and optical flux-limited sam- 
ple of AGN, and so in §4.4 we account for the spectro- 
scopic incompleteness to infer the AGN population to 
/o.5-iokeV < 1 X 10"^^ erg cm'^ and i^^ < 23. 

4.1. X-ray Flux Limit 

The first limit on the completeness is the target selec- 
tion, which is limited in both X-ray and optical fluxes. 
The initial selection includes all XMM targets with X- 
ray flux limits of 1 x 10~^^ erg cm~^ s~^ in 0.5-2 keV or 
6 X 10""'^^ erg cm~^ s~^ in the hard 2-10 keV band with 
optical counterparts of < 23.5. The X-ray flux limit 
means that we are complete in X-rays to all AGN with 
-^0.5-10 kcv > 3 X 10**^ erg (a classic AGN definition 
discussed in §3.2) at z < 1. 

4.2. Optical Flux Limit 

Our initial magnitude cut was i^g < 23.5, but this 
was designed to be more ambitious than the capabili- 
ties of Magellan/IMACS in 5-hour exposures. In Fig- 
ure [11] we show the spectral signal-to- noise (S /N) with 
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optical ij^^g magnitude for the targets from all IMACS 
exposures. The S/N was calculated by empirically mea- 
suring the noise in the central 6600-8200A region of each 
spectrum. The S/N generally correlates with the optical 
brightness, with some scatter attributable to varied con- 
ditions over the three years of observations. The outliers 
with high S/N and faint magnitude are all emission line 
sources where a strong emission line lies in the spectrum 
but outside the observed ijg filter range. The low-S/N 
and bright magnitude objects of the lower left may be 
highly variable sources or targets with photometry con- 
taminated by blending or nearby bright stars. The in- 
creasing number of unclassified targets (filled green cir- 
cles) in Figure[TT]shows that we do not identify all objects 
to iXb < 23.5. 

In Figure [T2I we show the completeness with i^g mag- 
nitude for the various classifications. We assume that the 
identified fractions have Poisson counting errors from the 
number of the given type and the total number of targets 
in each magnitude bin. The survey completeness to all 
targets remains at ^^90% to < 22. The identification 
fractions of emission line targets remains nearly flat a 
magnitude deeper than the absorption line galaxies, al- 
though the fractions of "bl" and "nl" objects decrease 
slightly from 22 < z^g < 23, within the noise. 

The completeness is not uniform for all types of 
objects: the fraction of identified broad and narrow 
emission line targets remains statistically constant un- 



for emission line spectra. We assume that the SPSS 
Type 1 composite spectrum (iVanden Berk et all 120011) 



and Type 2 composite spectrum (jZakamska et al 



2003) 



til 



AB 



22.5, while the fraction of absorption line tar- 



gets appears complete only to ijg ^ 22. Both narrow 
and broad emission lines generally exhibit two or more 
times the signal of their continuum, allowing for iden- 
tification even when the objects' broad-band magnitude 
and average S/N are low. Since different emission lines 
vary in strength, this also suggests that the identification 
of "bl" and "nl" may suffer from a redshift dependence 
(for instance, some redshifts may have only weak emis- 
sion lines in their wavelength range, while others include 
several strong lines). 

4.3. Redshift- Dependent Completeness 

The strongest redshift dependence in the spectra come 
from targets with only one strong emission line in their 
observed IMACS spectra. The "a" type objects are 
well-populated with absorption lines and have redshift- 
independent classifications, but emission line spectra 
may have only one line in the observed 5600-9200 A win- 
dow (see Figure ^ . The presence of only one emission 
line causes two problems: the redshift solution will be 
degenerate, and the line may be confused with noise if 
it is either narrow or broad but weak. The first prob- 
lem means that we can only assign Zconf = 2 and we 
may also assign the wrong redshift (bright targets are 
an exception, since a redshift can be assigned based on 
strong continuum features). The second problem means 
that we might completely miss the AGN designation and 
assign it a "?" classification. At faint S/N levels, the 
pattern of two emission lines is much easier to identify 
over the noise, and so the lower completeness to single- 
line objects may mean both lower redshift confidence and 
a lower identification threshold. 

We used Monte Carlo simulations to test the redshift 
and magnitude dependence of the survey's completeness 



each have infinite signal-to-noise, and degrade these spec- 
tra with Gaussian-distributed random noise to represent 
broad and narrow emission line spectra of varying ij^g 
magnitude. For each bin in magnitude, we calculate the 
median signal-to-noise of the observed spectra at that 
brightness, measured at both A < 8000 A and in the 
noisier region with sky lines at A > 8OOOA (our spec- 
tra typically have S/N about 16% worse at A > 8OOOA). 
Each artificial noise-added spectrum was then redshifted 
over several values and realized in the IMACS wavelength 
range (5600A-9200A). We then used the same idlspec2d 
redshift algorithm used on the data described in §3.1 to 
determine whether or not we would be able to assign the 
correct redshift with high-confidence (^conf — 3,4) for 
these artificial redshifted spectra (a redshift could not 
be determined if the emission lines were smeared out or 
if the spectrum could not be distinguished from noise or 
a different line at another redshift). We used 20 real- 
izations for each redshift and signal-to-noise bin. The 
fraction of artificial spectra with determined redshifts at 
a given redshift and signal-to-noise, with different seeds 
of randomly-added noise, forms an estimate of the com- 
pleteness. 

We found that the simulated completeness for narrow 
emission line spectra was 90% complete to «Jg ~ 23 
{S/N w 1.76 per pixel) for 2 < 1.3, with strong un- 
ambiguous lines (e.g., HaA6563, H/3A4861, [Oiii]A5007, 
[Oii]A3727). This is a magnitude fainter than the level 
of the average redshift completeness of the survey. At 
0.9 < z < 1.4, [Oii]A3727 is the only strong line, but it 
is bright enough that the redshift solution remains un- 
ambiguous even to ijg < 23. At z > 1.4 the [On] line 
shifts completely out of the wavelength range and no 
good emission lines remain. The additional blue Hec- 
tospec coverage is also useless at z > 1.4, since [On] re- 
mains redward of the upper 9200 A wavelength limit. We 
cannot identify narrow emission line ("nl") spectra at 
z > 1.4. 

The Type 1 AGN completeness has a more complex 
redshift dependence. As shown in Figure [6l in the red- 
shift ranges 0.4 < z < 1.9 and 2.3 < z < 2.9 only one line 
is present and the redshift solution may be degenerate. 
This is ameliorated by the ancillary MMT/Hectospec 
spectra which have broader wavelength coverage. Exam- 
ples of two objects with only one emission line in their 
IMACS spectra, but two emission lines in their Hectospec 
spectra, are shown in Figure [7] with accompanying dis- 
cussion in §3.3. Only 36% (104/288) of the broad emis- 
sion ("bl") spectra benefit from MMT/Hectospec cov- 
erage. We add this MMT/Hectospec corroboration to 
the unidentified targets in the simulations, and estimate 
the redshift completeness as shown in Figure [T31 We 
have lower redshift completeness in the redshift ranges 
0-5 ^ z < 1.5, and 2.3 ^ z < 2.6, where only one line 
is present (H/3, Mgn, or Cm]) and although we can re- 
liably classify as a broad line AGN ("bl") it is difficult 
to distinguish between the two redshift ranges. Without 
the degeneracies between redshift, the redshifts would be 
> 90% complete to S/N « 1.75 (per pixel) or i^g ~ 23). 

We do not test redshift dependence in identifying ab- 
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sorption line ("a") spectral types because these spectra 
are well-populated with absorption lines. At z > 1.3 the 
4000A break leaves the wavelength range, but otherwise 
each absorption line galaxy has the same aptitude for 
•^conf > 3 classification at z < 1.3. However, because 
the absorption line ( "a" ) spectra lack features which are 
of higher signal than their continua, we cannot identify 
them to the same low S/N levels as emission line spectra. 
So the incompleteness to absorption line ("a") spectra 
at z < 1.3 with 22 < z^g < 23 is not redshift depen- 
dent. Because we have high completeness to broad line 
AGN ("bl") and narrow emission line spectra ("nl") at 
this magnitude (excepting the redshift ranges described 
above), most of the unidentified targets at 22 < i^g < 23 
are probably absorption line galaxies. 

In summary, the sample has the following incomplete- 
ness outside of the flux limits: 

1. Type 1 AGN of 22 < < 23 at z - 0.8, z - 1.3, 
and 2 ~ 2.4 (completeness in these regions shown 
in Figure fTS]) . 

2. Type 2 AGN of all magnitudes at z > 1.4 

3. Absorption line galaxies of 22 < < 23 at z < 
1.3 (from §4.2), and absorption line galaxies of all 
magnitudes at z > 1.3 

We show the redshift distribution of all AGN (meeting 
one of the X-ray criteria in §3.2) in FigurelMl The uncor- 
rected redshift distribution is shown by the square sym- 
bols. We next attempt to describe the complete < 23 
flux-limited sample, correcting for the incompleteness of 
the four points above. 

4.4. Characterizing the Low- Confidence Targets 

We can only assign high-confidence (zconf = 3, 4) red- 
shifts for 72% (485/677) of the targets, leaving 88 spectra 
with low-confidence (zconf = 1,2 redshifts and 104 tar- 
gets of unknown spectral type (zconf = ~ljO). We char- 
acterize these 192 low-confidence and unclassified spec- 
tra using the photom etric classifications and redshifts of 
ISalvato et al.l (|2008f l. which take ad vantage of the exten- 
sive photometry of COSMOS fCap ak et al. l [200l . The 
photometric redshift algorithm finds a best-fit redshift 
and classification by matching a set of 30 templates to the 
IR (IRAC), optical (Subaru), and UV (GALEX) photo- 
metric data of e ach object. Th e templates are described 
in full detail in ISalvato et al.l (j200S) and are available 
upon request^^. The photometric redshift technique was 
calibrated upon the spectroscopic redshifts we present 
for the 485 spectra of high redshift confidence, and has 
a precision of (Taz/(1 + z) < 0.02 with < 5% of targets 
as significant outliers at z < 4.5. 

The photometric redshift templates rely on multiwave- 
length fitting from IR to UV wavelengths, and so the 
photometric classifications can separate AGN-dominated 
(which we designated "unobscured" ) and host-dominated 
(which we designate "obscured") AGN types. However, 
the photometric classifications do not distinguish well be- 
tween our absorption line spectra ("a" types) and nar- 
row emission line spectra ("nl" types), although they can 
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separate unobscured broad line AGN ("bl" types) from 
obscured AGN ("a" and "nl" types). We must assume 
population fractions of absorption line and narrow emis- 
sion line spectra from the photometrically classified ob- 
scured objects using the known fractions from the high- 
confidence spectroscopy. In Figure [121 the fraction of 
narrow emission ( "nl" ) spectra does not decrease appre- 
ciably to jjg ~ 22.5 and almost all of the unknown ob- 
jects can be assumed to be absorption line ("a") types. 
We also know from §4.3 above that we are incomplete to 
"nl" spectra at z > 1.4, and the spectroscopically unclas- 
sified targets at z > 1.4 probably follow the ~2:1 ratio 
of narrow emission ( "nl" ) to absorption ( "a" ) types we 
find at lower redshifts in §3.1. So we assume that all 
photometrically classified unobscured AGN correspond 
to our broad emission ("bl") type, and assume fractions 
of absorption ( "a" ) and narrow emission spectra ( "nl" ) 
as follows: (1) for z < 1.4, all are "a" types, and (2) for 
z > 1.4, 2/3 are "nl" types and the remainder are "a" 
types. 

Figure [15] shows the photometric redshift distribution 
for the 192 low confidence and unclassified spectra. Most 
of the photometric redshifts fall into one of the three in- 
completeness categories shown above in §4.3. We will use 
this redshift distribution to characterize the demograph- 
ics of the complete flux-limited sample in §5. 

We can also use the absolute magnitude distribution of 
the targets with secure spectroscopic redshifts in Figure 
1161 to make a qualitative assessment of the unidentified 
targets. The dashed lines mark i^g = 22 and ijg — 23. 
We will assume that the absolute magnitude distribution 
for narrow emission ("nl") and absorption ("a") objects, 
which peaks at Mi ~ —22, does not change with red- 
shift. Objects of this absolute magnitude distribution 
should be detected to z ~ 2, but there are no narrow 
emission ("nl") spectra detected at z > 1.4, and no ab- 
sorption ("a") spectra detected at z > 1.3. So many 
of the unidentified targets are probably z > 1.4 "nl" and 
z > 1.3 "a" type objects. The bright tail of the Mi distri- 
bution for "a" and "nl" types at z > 1.3/z > 1.4 also has 
i^g < 22, and so these missing z > 1.3/z > 1.4 targets 
may account for the unidentified targets at i^g < 22. 
In addition, most obscured AGN have Mi ~ —22, which 
hes within 22 < i^g < 23 at z > 1.3/z > 1.4, suggest- 
ing that these objects may be most of the unidentified 
22 < i^g < 23 targets. This qualitative assessment con- 
firms the characterization of the unknown spectral types 
using photometric redshifts. 

5. DISCUSSION 

5.1. Demographics 

Figure [HI indicates that the X-ray flux limit includes all 
-^0.5-iokcV > 3 X lO'^^ erg s^^ AGN at z < 1. This means 
we are nearly complete to all X-ray AGN (as defined in 
§3.2) at z < 1, since almost all objects that meet one of 
the X-ray criteria also meet both. We can additionally 
see in Figure [16] that we observe all but the faint tail of 
the Mi distributions of obscured and unobscured AGN 
types to z < 1, so long as we use the simple corrections 
of §4.4 to characterize the sample to i^g < 23. This 
allows us to characterize the complete z < 1 X-ray AGN 
population. 

In Figure [14] we show the number of each AGN type 
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with redshift. This includes only the 432 high-confidence 
X-ray AGN as defined by the X-ray criteria. We find raw 
fractions of broad emission line (56%), narrow emission 
line (32%), and absorption line (12%) over all redshifts, 
which rough ly agree wi th other wide-area X-ray surveys 
(iFiore et all I2b03- Silverman et all l2005l : lEckart et all 
120061 iTrump et al. I l2007fl . To characterize the complete 
i^B < 23 and /0.5-iokoV > 1 x 10~^^ erg s~^ cm~^ sam- 
ple, we include the 106 ijg < 23 targets with bad spec- 
troscopy and photometric redshifts that satisfy the X-ray 
AGN criteria. The corrected fractions of i\g < 23 tar- 
gets at all redshifts include 57% broad emission line, 25% 
narrow emission line, and 18% absorption line AGN. 

5.2. Obscured to Unobscured AGN Ratio 

The ratio of obscured to unobscured AGN can help 
determine the properties of the obscuration which hides 
nucl ear activit y. I n the simplest unifica tion mod- 
els (|Antonuccj I1993I : lUrrv fc Padovanil 1 1995( 1 . obscura- 
tion depends only on the orientation and should re- 
main independent of luminosity and redshift. How- 
ever, since we know that galaxies at higher redshifts 
have more dust than local galaxies, then one might 
expect the ratio of obscured to unobscured AGN to 
depend on redshift if AGN host galaxy dust plays a 
role i n obscuration (e.g.. iBallantvne. Everett, fc Murravl 
I2OO6D . And if the obscuring dust (or its sublimation ra- 
dius) is blown out further by more luminous accretion 
disks ([Lawrence fc Elvis! IT981 ILawrencel[T99lt ISimpsonI 
[2OO5I) . then one might expect the ratio to decrease 
with increasing luminosity. Some models of the X-ray 
background prefer ratios which suit these physical de- 
scriptions, predicting an increasing ratio of obscured 
to unobscured with in creasing redshift and d e creas- 
ing luminosity (IBallant vne. Everett, fc Murravl 120061 : 
iTreister fc Urrvl l2006l ). Deep X-ray obse rvations con- 
firm t hat the ratio depeiids on luminositv (ISteffen et al. 
20041 iBarger et af] l2005t ITreister. Krolik fc DuUemond 
2008() . Some observations additionally suggest redshift 



iHasingeil [200S1 , but other authors claim that redshift 
evolution is neither nece ssary in the models nor signifi- 
cant i n the observatio ns (jUeda et al.l 120031 : lAkvlas et al.l 
I2OO6I: iGilh erani200l . 

We derive the obscured to unobscured AGN ratio with 
redshift in Figure [171 Here "obscured AGN" refers to 
both narrow emission line ("nl") and absorption line 
( "a" ) AGN meeting the X-ray criteria of §3.2, and "unob- 
scured AGN" includes all broad-hne ("bl") AGN. To the 
limit of the survey at z < 1, our average ratio is 3:1 ob- 
scured to unobscured AGN. We additionally separate the 
AGN into X-ray luminous and X-ray faint (in relation to 
the median X-ray luminosity, £mcd — 1-32 x 10*^ cgs) in 
the bottom panel of Figure [TTl The ratio of obscured to 
unobscured X-ray faint AGN appears to be much higher 
than the ratio for X-ray bright AGN, and additionally 
seems to increase with redshift. 

We test the ratio for dependence on redshift and lumi- 
nosity using logistic regression a useful method for deter- 
mining how classification depends upon a set of variables. 
It is commonly used in biostatistical applications, where 
one expects a binary response (for instance, a patient 
might live or die) based on a set of variables. Logistic re- 



gression considers each data element independently, and 
is therefore more effective than significance tests which 
bin the data. An ex cellent review of logistic regression is 
found in|^ (IT99I. We use the method here to learn if 
the likelihood for an AGN to be classified as obscured or 
unobscured (a binary response) depends on observed X- 
ray luminosity and/or redshift. Logistic regression solves 
for the "logit" (the natural logarithm of the odds ratio) 
in terms of the variables as follows: 



In 



Pr{G = l|L,z) 
Pr{G = 2\L,z) 



(3o + Pi ln(l + z)+l32 HLx/l^n 

(3) 

Here G = 1 means an AGN is classified obscured, and 
G — 2 means an AGN is classified unobscured. We use 
log(l -I- z) and \og{Lx /lO'^^) as the dependent variables 
instead of z and Lx for numerical stability. Then the 
logit, as the logarithm of the ratio of the probabilities, is 
just the logarithm of the obscured to unobscured ratio. 
We solve for the coefficients using the Newton-Raphson 
method and estimate errors by bootstrapping, calculat- 
ing the standard deviation on the coefficients with 1000 
random subsets of the true data. We find the coefficients 
to be 



l3o = 1.73±0.67; /3i = 3.83±1.57; P2 = -0.69±0.19 (4) 

In other words, the obscured/unobscured AGN ratio 
increases with redshift at 2.4cr significance and decreases 
with observed X-ray luminosity at 3.6cr significance. We 
can write the power-law equation of the expected ratio 
for a given luminosity and redshift as: 



Pr (obscured) 
Pr(unobscured) 



oc 5.6(1 + z)3-8(ix/10^2)- 



(5) 



The curves from this logistic regression model are 
shown in Figure [T71 In the bottom panel, the red line 
shows the power-law relation (Equation 5) computed us- 
ing Lx = iiow, where Liow is the median luminosity 
from only those AGN with _Lo.5-io kcV < -^mcd- Sim- 
ilarly the blue line represents the relation for higher- 



luminosity AGN of L0.5- 



10 kcV 



Note that the 



data as binned in Figure [T7| shows less signal than the 
independent data used in the logistic regression fit, and 
so the fit of the power-laws shown should not be judged 
by the basis of their by-eye match to the binned data. 
It is worth noting, however, that high-luminosity sources 
seem to evolve much more weakly with redshift. This is a 
natural consequence of the power-law nature of Equation 
5: when Lx is large, the obscured/unobscured ratio be- 
comes small, and so it appears only to weakly evolve with 
redshift on a linear scale. Our dependence of obscuration 
on redshift and lumino sity are bot h consistent with re- 
cent w ork by both Treis ter fc Urrvl ({2006) and Hasingl 
(|2008l ). with the obscured fraction about 4 times higher 
at low luminosity than at high luminosity and about 2 
times higher at z 1 than at z ~ 0. 

The trend with observed X-ray luminosity can be ex- 
plained in several ways. It may be that obscured AGN 
are simply more absorbed in the X-rays, such that their 
intrinsic X-ray luminosities are significantly higher than 
their observed. Then the apparent lack of obscured AGN 
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at higher X-ray luminosities might be only an observed 
effect, and not an intrinsic physical effect. But if the 
intrinsic and observed X-ray luminosities are not signif- 
icantly different in these obscured AGN, then the lumi- 
nosity dependence indicates that more luminous AGN 
have less obscuring material. The luminosity may de- 
crease the opening angle of obscuration by causing dust 
sublimation to occur at larger radii. 

The presence of more obscured AGN at z ~ 1 is not 
likely to indicate physical evolution in AGN, since AGN 
at similar luminosities at z ~ 1 and z ~ are not ob- 
served to have different ph ysical properties li ke black hole 
mass and accretion rate (iKellv et al.l l2008^ or spectral 
energy distribution (IVignali et al.l 20031: Richards et al.l 



|2006[ iHopkins. Richards, fc Hernauistll2007f ). However, 
galaxies at z ~ 1 show significantly more star formation, 
gas, and dust than galaxies at z ~ 0, and so the increase 
of obscuration with redshift may be explained by host 
gas/dust o bscuration o f the A GN central engine. Indeed, 
models by iBallantvn^ (120081 ) show that star formation 
can effectively obscure AGN while producing both the 
observed luminosity and redshift dependen ce of the ob- 
scured/unobscured ratio. iBallantvng (|2008f l additionally 
show that starburst-driven obscuration should be easily 
distinguished from AGN-heated dust by future Herschel 
100 fim surveys. 

It is important to note that our definition of "ob- 
scured" includes only moderately X-ray obscured AGN. 
We are not sensitive to Compton-thick and other heavily 
X-ray obscured AGN, and so may be significantly under- 
estimating the obscure d AGN population (Dad di et al.l 
I2007t iFiore et al.ll2008[ ). Logistic regression reveals sta- 
tistically significant evidence for redshift evolution and 
dependence on X-ray luminosity in the optically ob- 
scured/unobscured ratio, but mid-IR surveys may re- 
veal different dependencies by including heavily obscured 
AGN missed in X-rays. 

6. CONCLUSIONS AND FUTURE PROJECTS 

We present optical spectroscopy for 677 X-ray targets 
from COSMOS, with spectra from Magellan/IMACS, 
MMT/Hectospec, and archival SDSS data. The spec- 
troscopy is uniformly complete to < 22. By using 
photometric redshifts for the bad spectra, we addition- 
ally characterize the sample to < 23, and we show 
that this optical limit, along with our X-ray flux limit, 
allows us to characterize a solely volume-limited sample 
of all (obscured and unobscured) X-ray AGN at z < 1. 
We provide evidence that at z < 1, the ratio of obscured 
to unobscured AGN increases with redshift and decreases 
with luminosity, where the redshift dependence is of mod- 
erate statistical significance (2.4(7) and the luminosity 
dependence is of higher statistical significance (3.6(t). 

Despite such leverage in the sample presented here, the 
observations of AGN in COSMOS are by no means com- 
plete. We were only able to target 52% of the available 
< 23.5 XMM targets, and we hope to include the re- 
mainder of targets in future spectroscopic observations. 
Some of these targets were observed on Magellan/IMACS 



and MMT/Hectospec in March 2008, and many of the 
other XMM targets without spectra will be observed with 
VLT/VIMOS (at 5600-9400 A) as part of the zCOS- 
MOS galaxy redshift survey (jLillv et al.l l2008f) . The 
zCOSMOS survey will additionally target > 23.5 
XMM targets which are too faint for Magellan/IMACS. 
It is also possible to study fainter X-ray sources, since 
the 0.5-2 keV Chandra observations in COSMOS go to 
2 X 10~^^ erg cm~^ s~^ in the central 0.8 dcg^, five times 
fainter than the XMM observations used here. Optical 
identification of these sources are still ongoing, but the 
Chandra data are expected to reveal twice as many X- 
ray targets as the XMM-selected targets presented here. 
We will additionally use the previously observed spectra 
of radio and infrared selected AGN candidates to study 
Compton-thick and other X-ray faint AGN, 

Future work will also use the bolometric studies made 
possible by the deep multiwavelength coverage of COS- 
MOS. We plan to further study the evolution of obscura- 
tion with more fundamental physical quanti ties like bolo- 
metr ic luminosity. A companion paper (jTrump et al.l 
l2008l ) presents virial black hole mass estimates for the 
Type 1 AGN presented here and suggests that it is dif- 
ficult to form a broad line region belo w a critical accre- 
tion rate, as sug gested pre viously by N icastro &: Elvid 
(|200g) and iKoUmeier et al.l (|2006). This concept, com- 
bined with the luminosity evolution of the obscuration 
presented here, suggests that models of the AGN central 
engine must include a prescription where the amount of 
obscuring material decreases with increasing luminosity, 
accretion rate, or both. 
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TABLE 1 

COSMOS Observation Log of the IMACS Spectroscopic Observations 



IMACS 


Center (J2000) 


Observation 


Exposure 


Number of 


Field 


RA (hh:mm:ss) 


Dec (dd:mm:ss) 


Year 


(hours) 


Spectra Extracted 




09:58:24 


02:42:34 




4 52 


oo 


2 


09:59:48 


02:42:30 


2007 


4.00 


46 


3 


10:01:06 


02:42:38 


2006 2007^ 


6.00 


75 


4 


10:02:33 


02:42:34 


2006 


5.33 


37 


5 


09:58:26 


02:21:25 


2006,2007=^ 


6.00 


56 


6 


09:59:47 


02:21:25 


2005 


6.90 


43 


7 


10:01:10 


02:21:25 


2005 


6.03 


48 


8 


10:02:36 


02:21:29 


2006 


5.10 


42 


9 


09:58:25 


02:00:13 


2006 


5.03 


27 


10 


09:59:47 


02:00:17 


2005 


6.64 


35 


11 


10:01:10 


02:00:17 


2005 


4.67 


39 


12 


10:02:37 


02:02:05 


2005 


4.77 


37 


13 


09:58:24 


01:39:08 


2006 


2.67 


7 


14 


09:59:47 


01:39:08 


2006 


5.33 


23 


15 


10:01:10 


01:39:08 


2005 


3.63 


25 


16 


10:02:33 


01:39:08 


2005 


3.93 


28 



Fields 3 and 5 were observed for one hour in 2006 and five hours in 2007, for six total hours of exposure. 



TABLE 2 

COSMOS XMM Optical Spectroscopy Catalog 



Object Name 


RA {J2000) 

cict^rccs 


Dec (J2000) 

(J2000)'' 


i+ 

'CFHT 

AB mag 


S/N 


texp 
sec 


Type 


z 


(7^ 


b 

"^conf 


SDSS J095728.34+022542.2 


149.3680700 


2.4283800 


19.64 


7.00 





bl 


1.5356 


0.0015 


4 


COSMOS ,1095740.78+020207.9 


149.4199229 


2.0355304 


21.55 


17.92 


19200 


bl 


1.4800 


0.0028 


4 


SDSS ,1095743.33+024823.8 


149.4305400 


2.8066200 


20.43 


3.37 





bl 


1.3588 


0.0020 


4 


COSMOS J095743. 85+022239.1 


149.4327000 


2.3775230 


23.40 


1.33 


18000 


nl 


1.0192 


0.0002 


1 


COSMOS J095743. 95+015825.6 


149.4331452 


1.9737751 


21.91 


1.54 


19200 


a 


0.4856 


0.0030 


1 


COSMOS J095746. 71+020711. 8 


149.4446179 


2.1199407 


20.78 


11.37 


19200 


bl 


0.9855 


0.0002 


4 


COSMOS J095749. 02+015310.1 


149.45426.38 


1.8861407 


20.36 


13.68 


19200 


Ilia 


0.3187 


0.0002 


4 


COSMOS J095751.08+022124.6 


149.4628491 


2.3568402 


20.73 


3.91 


18000 


bnl 


1.8446 


0.0001 


4 


COSMOS J095752. 17+015120.1 


149.4673623 


1.8555716 


21.08 


7.31 


19200 


bl 


4.1744 


0.0005 


4 


COSMOS J095753.44+024114.2 


149.4726733 


2.6872864 


22.18 


0.95 


11160 


bl 


2.3100 


-1.0000 


1 


COSMOS J095753.49+024736.1 


149.4728835 


2.7933716 


21.96 


4.76 


11160 


bl 


3.6095 


0.0128 


4 


SDSS J095754. 11+025508.4 


149.4754500 


2.9189900 


19.45 


6.09 





bl 


1.5688 


0.0022 


4 


SDSS J095754.70+023832.9 


149.4779200 


2.6424700 


19.35 


8.04 





bl 


1.6004 


0.0015 


4 


SDSS J095755.08+024806.6 


149.4795000 


2.8018400 


19.41 


8.66 





bl 


1.1108 


0.0017 


4 


COSMOS ,1095755.48+022401.1 


149.4811514 


2.4003076 


21.26 


19.89 


18000 


bl 


3.1033 


0.0003 


4 


COSMOS ,1095756.77+024840.9 


149.4865392 


2.8113728 


20.81 


11.86 


11160 


bl 


1.6133 


0.0098 


3 


COSMOS J095757. 50+023920.1 


149.4895683 


2.6555795 


20.30 


11.17 


11160 


nl 


0.4674 


0.0002 


2 


SDSS J095759.50+020436.1 


149.4979100 


2.0766900 


18.98 


14.57 





bl 


2.0302 


0.0016 


4 


COSMOS J095800. 41+022452.5 


149.5017000 


2.4145710 


22.57 


3.53 


18000 


bnl 


1.4055 


0.0001 


4 


COSMOS J095801. 34+024327.9 


149.5055777 


2.7244216 


20.66 


9.67 


11160 


nla 


0.3950 


0.0010 


1 


COSMOS J095801. 45+014832.9 


149.5060326 


1.8091427 


21.96 


1.79 


9600 


bl 


2.3995 


0.0002 


4 


COSMOS J095801. 61+020428.9 


149.5067217 


2.0746879 


22.18 


5.46 


19200 


bl 


1.2260 


-1.0000 


1 


COSMOS J095801. 78+023726.2 


149.5074058 


2.6239318 


17.79 


52.98 


11160 


star 


0.0000 


0.0000 


4 


COSMOS J095802.10+021541.0 


149.5087524 


2.2613900 


21.01 


3.75 


3600 


a 


0.9431 


0.0050 


3 



The RA and Dec refer to the optical counterpart of the X-ray source, which is where the slit was centered.'' From empirical measurements, the redshift confidence 
was found to correspond to correct redshift likelihoods of 97%, 90%, 75%, and 33% for Zconf = 4, 3, 2, 1, respectively. The redshift confidences are fully explained 
in §3.1. 
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TABLE 3 
Targeting and Redshift Yields 





16-Field 




Per Field 




X-ray Sources 


Total 


Minimum 


Maximum 


Median 


All Sources 


1640 


68 


145 


105 


i+B < 23.5 


1310 


55 


110 


86 


Targeted 


677 


9 


74 


38 


Classified (^^onf > 0) 


573 


9 


63 


30 


^conf = 3,4 Redshifts 


485 


8 


53 


26 


^conf = 3,4 with Hcctospec 


117 


1 


27 


6 


Zconf =3,4 with SDSS 


76 


2 


12 


4 
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Fig. 1. — The X-ray depth and survey size of various deep X-ray AGN surveys, along with the X-ray and optical flux for targets 
in COSMOS. At left, symbol sizes indicate each survey's number of X-ray point sourc es: open indicate all source s, and filled indi- 
cate those with optical spectroscopy. References for the sur veys are as follows: AEGIS f Davis et al."2007'), AGES (Brand et al.''2006), 
CDF-N (Alexander et al. 2003; Bargcr ct al. 2003), CDF-S (ILuo et al. 2008), CHAMP (Kim et al. 2004; Green et al. 2004), CLASXS 
(lYang -t al. 2004), HELLAS2XMM (Fiore ct al. 2003: Co cchia et al.112007 '1. LALA (Wane et al. 2004), RDS/XMM (Lehmann et al. 2001), 
RMS (Hasinger et al. 2005), and SXDS (Ucda ct al. 2008). At right, the crosses represent all XMM point sources from Brusa et al. (200^) 
and the gray shaded area represents the "AGN locus" of — 1 < log(/x //o) < 1 iMaccacaro et al... 1988 ). The COSMOS Chandra data (not 
presented here) go four times deeper in the central 0.8 deg^ , doubling the number of COSMOS point sources. 
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Fig. 2. — Maps of X-ray sensitivity for each of the three years of IMACS observing. The top left panel shows the XMM depth and 
IMACS pointings for the first year, the top right shows the second year cumulative depth and pointings, and the lower left shows the third 
year cumulative depth and pointings. Since the XMM observations were ongoing during the spectroscopy campaign, we chose each year's 
IMACS pointings from the regions of greatest XMM uniformity and depth lacking previous spectroscopic observations. 
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Fig. 3. — Spectroscopic observations of the 2 deg'^ COSMOS area. X-ray targets with IMACS spectra are shown as crosses, those with 
MMT spectra are diamonds, and those with SDSS spectra are squares. The 16 tiled IMACS pointings are shown as boxes of 22'30" X 21' 10" 
and are shaded according to their exposure time. The two 1-deg diameter MMT pointings are shown as circles. COSMOS also includes 
deeper Chandra coverage, not used here, over the field's central square degree (fields 6, 7, 10, and 11, with portions of the other 8 fields). 
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Fig. 4. — Five examples of IMACS spectra with broad emission lines. The dominant line species are labeled in each spectrum and bad 
pixels are omitted. The first two objects are Type 1 AGN ("bl") with tiic iiigiicst rcdsliift confidences, the third is a high-confidence AGN 
with both narrow and broad emission ("bnl"), and the bottom two are Type 1 AGN ("bl") with uncertain redshifts. We discuss these 
objects in §2.2. 
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Fig. 5. — Five more examples of IMACS spectra, including four targets with narrow emission lines and one absorption line galaxy. 
The prominent absorption and emission features are labeled. The third target is a hybrid "nla" object with both narrow emission and 
absorption lines. The first four objects have the highest rcdshift confidence, while the bottom target has an extremely uncertain redshift, 
calculated from a single emission line which may be solely due to noise. The first is a starburst galaxy by its X-ray emission, while the 
second, third, and fourth spectra are all AGN which meet both of the X-ray criteria of §3.2. We discuss these objects in §2.2. 
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Fig. 6. — The observed wavelengths of prominent broad emission Unes with redshift. The spectral ranges of MMT/Hectospec and 
Magellan/IMACS are shown at the top. The broad emission lines observed at a given redshift can be found by drawing a horizontal line 
between the wavelength limits: the solid lines of broad emission peak intersecting that redshift line would be present in the spectrum. 
The narrow wavelength coverage of IMACS means that only one broad line is present in the shaded redshift ranges 0.4 < z < 1.9 and 
2.3 < z < 2.9, so that spectra with low S/N may be assigned z^onf = 2 because they have degenerate redshift solutions. The extended 
wavelength coverage of Hectospec allows us to resolve the degeneracies and assign Zconf = 4. 



COSMOS AGN 



21 



3.0 
2.5 
2.0 
1.5 
1.0 

^ 0.5 

•.^ 0.0 
3.0 

u 

> 2.5 

^ 2.0 

."^ 1-5 

C 0.5 

; 0, 



095801.45 + 014832.9 bl i^g =21.95 z = 2.40 z,„„,=4 



CHI] 



IMACS 



_Lya 



NV/Sill 



CIV 




MMT 





100149.00 + 024821.8 bl =22.99 z=1.61 z„„,=4 IMACS 

Mgll 




4000 



5000 



6000 7000 
Wavelength (A) 



8000 



9000 



Fig. 7. — Two X-ray targets with both IMACS and Hectospec spectra. In both cases, the IMACS wavelength range only includes one 
broad emission line and so has a degenerate redshift solution. The additional blue Hectospec coverage resolves the degeneracy and allows 
us to assign these objects z^onf = 4- The fourth panel shows that our Hectospec flux calibration can cause errors in spectral shape at red 
wavelengths, although this does not ailect our redshift solutions. 
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Fig. 8. — The three templates used in the c lassification and rcdshift determination scheme. The broad line AGN template is the SDSS 
quasa r composite of [Va ndcn Be rk et al.l 1 1200 II ). the narrow emission line template is the SDSS Ty pe II AGN composite o f Zakamska ct al] 
II2003I ) . and the absorption line red galaxy template is the composite of the SDSS red galaxy sample IIEisenstein et al.l2001l ). The wavelength 



coverages of the templates were sufficient for the entire redshift range (and the corresponding observed wavelength ranges) of the sample. 
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Fig. 9. — The observed 0.5-10 keV X-ray luminosities for z^onl > 3 objects with redshift. We label objects classified as "bl" (Type 1 
AGN) with black crosses, "nl" with blue diamonds, and "a" with red squares. The dashed line shows the survey's limiting luminosity from 
the XMM flux limit. The AGN luminosity cut of io.5— lOkcV = 3 X 10"*^ is drawn as a solid line: all "nl" and "a" above this line are AGN. 
Objects below this line, however, are not necessarily inactive: two Type 1 AGN are less luminous, and the luminosity limit is conservative. 
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Fig. 10. — The X-ray flux vs. the ijg magnitude for the z^onf ^ 3 objects. The AGN locus of —1 < log fx / fo ^ 1 ("Maccacar o et aTl 
119881 . , see also Equation 2) is shown by the gray shaded region, along with approximate boundary between quiescent and star-forming 
galaxies at log fx / fo = ~2 IjBauer et al.|[20 04y Black crosses are targets classified "bl" (Type 1 AGN), blue diamonds are "nl," and red 
squares are "a" objects. We additionally mark all targets of Lo.5— lOkcV > 3 X lO**^ with black x's. We consider targets either in the AGN 
locus or with Lx > 3 X 10*^ to be AGN: this includes all of the "bl" spectra and all but 53 of the "nl" and "a" spectra. 
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Fig. 11. — The signal-to-noise (S/N) and optical magnitudes for the X-ray targets. Crosses indicate broad emission line spectra, 
diamonds are narrow emission line spectra, squares are absorption line spectra, and filled circles are unclassified objects. The S/N and 
optical magnitude are correlated, with scatter from varying conditions over three years of observing. The number of unidentified objects 
increases greatly at > 22, although we still identify emission line spectra at the faintest magnitudes. 
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Fig. 12. — Cumulative completeness by classification type versus optical ij^g magnitude. Each region (shaded or unshaded) indicates the 
relative fraction of broad emission line objects ("bl"), narrow emission line objects ("nl"), and absorption line ("a") spectra. Targets with 
too low S/N to venture a classification axe represented in the upper "?" region. Error bars on the points above each region are calculated 
assuming that, in each magnitude bin, both the number of each class and the total number have associated Poisson counting errors. The 
total completeness for each classification is ~90% to ijg < 22, although we can correct for the incompleteness of each spectral type to 

4b < 23. 
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Fig. 13. — Our estimated completeness to assigning 2fconf > 3 redshifts for broad emission line spectra (Type 1 AGN). We used Monte 
Carlo simulations with 20 different spectra with Gaussian-distributed noise for each of 4 values of S/N and 20 redshift bins. Each point 
represents the number of simulated spectra assigned Zconf ^ 3, with an additional 36% of the bad Zconf < 3 simulated spectra based on 
the partial MMT/Hectospec coverage (since 104/288 observed Type 1 AGN had supplemental Hectospec spectra). Each signal-to- noise is 
associated with an it-n magnitude according to the median values in Figure 1111 The redshift ranges of lowest completeness correspond 



to observed wavelength ranges with only one emission line, as detailed in 
complete are identified as "bl" objects but with degenerate spectra. 



j4.3. Almost all of the simulated spectra to which we are not 
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Fig. 14. — The redshift distributions of broad line ("bl") AGN, narrow emission line ("nl") spectra, and absorption line ("a") galax ies. 
The raw distributions are shown with squares, while the distributions adjusted for the incompleteness (see §4.3 & 4.4 and Figure [T5)l are 
shown with diamonds. We do not correct the "nl" and "a" types where there is no spectroscopic data at z > 1.4. 
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Fig. 15. — The photometric redshifts for the spectra without high-confidence redshifts. Square symbols show all 192 objects, and 
diamonds show the 146 ijg < 23 objects. The spectral type for these objects comes from the template used for the photometric redshift, 
with "a" and "nl" fractions estimated as described in §4.4. We use the photometric redshifts and classifications to characterize the complete 
i < 23 X-ray AGN sample. 
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Fig. 16. — The absolute magnitude of the AGN (meeting the X-ray criteria of —1 < log fx / fo < 1 or io.S— lOkeV > 3 X 10*^ ) with 
redshift. The solid line indicates an arbitrary quasar/Seyfert boundary at Mi = —23, and the dashed lines indicate fluxes of ijg = 22 and 
= 23. Symbols for object types are as in Figures llll and [T2l with crosses for "bl" or "bnl" objects (Type 1 AGN), diamonds for "nl" 
and "nla" objects (Type 2 AGN), and squares for "a" objects (optically obscured AGN). 
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Fig. 17. — The ratio of obscured to unobscured AGN with redshift. We define obscured AGN as spectra with narrow emission or absorption 
lines ("nl," "a,", and "nla") that meet the X-ray AGN criteria of §3.2, while unobscured AGN are all broad-line ("bl" and "bnl") spectra. 
In the top panel the raw fractions are shown in gray, while the corrected fractions based on the incompleteness (characterized in §4.3 
& 4.4) are shown by the black diamonds. The bottom panel shows the ratios for AGN fainter and brighter than the median luminosity 
^med = 1-32 X 10*^ cgs, with the Lq.s— lO keV < ^med ratio as squares, and the Lq.s-io koV > ^med ratio as triangles. The errors associated 
with each point assume that the numbers of objects observed in each redshift bin are Poissonian. Logistic regression analysis shows that 
the ratio of obscured to unobscured increases with redshift and decreases with luminosity, as shown by the best-fit power-laws. 



